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This paper presents a graphical method of determining the groundwave attenuation 
over a spherical homogeneous earth in the far diffraction region. The curves are applicable 
to either vertical or horizontal polarization and to any combination of effective earth's 
radius, electromagnetic ground constants, frequency, path distance, and antenna heights. 
A criterion is given that indicates the met hod may be used not only for far beyond line-of- 
sight paths but, in many practical situations, at line-of-sight or even slightly within. 
Examples illustrating the use of the formulas and curves are included. 



1. Introduction 

In a previous paper [Vogler, 1961] a simplified 
graphical method was presented for calculating smooth 
homogeneous earth groundwave attenuation in the 
far diffraction region assuming horizontal polarization 
of the radio waves. In the present paper this met hod 
will be extended to include the case of vertical 
polarization. 

The far diffraction region is defined here as that 
region in which the diffracted field intensity may be 
determined by the first term of the Van der Pol- 
Bremmer residue series [Bremmer, 1949]. This 
region extends from near the radio horizon to all 
greater distances. In some situations the first 
term provides a valid approximation to the diffracted 
field even at points slightly within line-of-sight, and 
a criterion for determining the minimum distance 
is given in section 2. 

Since this paper is meant as a practical aid in 
obtaining diffraction calculations, no mathematical 
derivation of the functions discussed later is included. 
Detailed explanations of the residue series are given 
in the work by Bremmer [1949] and also in the 
introduction to the CCIR Atlas [Atlas, 1955]. As 
far as possible, the notation used is that established 
by Norton [1941] when identical quantities are 
considered; e.g., the parameters K and b°. The 
three functions C\, F, and 6? from which the attenu- 
ation is calculated were defined in such a manner as 
to facilitate graphical interpolation for any combi- 
nation of frequency and distance. Numerical evalu- 
ations of these functions were obtained through the 
use of an electronic computer. 

2. Groundwave Attenuation 

The attenuation A t (measured in decibels) relative 
to a free space inverse distance field (see, for example, 



(18)), may be divided into four terms containing, 
essentially, the distance dependence, two antenna 
height dependences, and a term depending on the 
electric constants of the ground: 



where 



A t =G(x )-F(x 1 )-F(x 2 )-C l (K } b°), (1) 



^i, 2 =ftC?/Uu, Co=(4/3*) J , (2) 



and /mhz denotes the radiofrequency in megahertz. 
da is the total arc distance between antennas, and 
rfi and d 2 are distances from the transmitting and 
receiving antennas to their radio horizons (see 
fig. 1). For antenna heights hi and h 2 not too large, 
the horizon distances may be approximated for a 
linear atmosphere with no ducts by 



d U2 ^2(ka)h lt2 
where ka is the effective radius of the sphere. 



(3) 




Figure 1. The geometry for smooth earth diffraction. 
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In practical applications the factor k is used to 
account approximately for the effects of atmospheric 
refraction [Burrows and Atwood, 1949] and is defined 
in terms of the gradient with respect to height, 
dn/dh, of the atmosphere's refractive index n: 



h 



L \n dh/ h =o. 



J' 



(4) 



It is usual to assume a standard "4/3 earth" 
atmosphere in presenting diffraction curves [Atlas, 
1955; Wait and Howe, 1956], and this procedure is 
followed here. Thus, the C appearing in (2) is 
merely a normalizing factor which for the standard 
atmosphere (fc=4/3) is equal to unity. Choosing a 
reference radius a=6373 km, (3) may now be 
expressed as 



di, 2 (km)~4.1225VCo" 3 Ai.2(m), 



(5) 



where h is in meters. The quantities x , x h and x 2 
defined by (2) can be considered as parametric 
distances; x referring to the total separation between 
antennas, and x u x 2 referring to the two antenna 
radio horizon distances. 

The parameters j3 , K, and b° in (1) and (2) 
depend on the polarization of the wave and the 
relative dielectric constant e and conductivity a 
of the ground. Assuming a homogeneous ground in 
which e and a are everywhere constant, K and b° 
are denned as follows for the two cases of vertical 



and horizontal polarization. Vertical polarization: 
K 9 =1.7778X10-*CJ££ [-^==J 2j ^ 

b° v =2 tan" 1 (jVtan" 1 (^J> (6b) 

horizontal polarization: 

K h =1.7778X10- 2 C Ml[(e-l) 2 +sT u \ (7a) 

b° h ^lSC^-tan- 1 (^~\ (7b) 

where s depends on the ground conductivity in 
mhos per meter and radiofrequency, and is defined 
as 

s=1.8X10 4 (r(mhos/m)//MHz. (8) 

Note that for a(mhos/m)// M Hz>>(e/2) X10" 4 : 

X f -2.385CiV5 : //!&, «~0, (9) 

^-1.325X10- 4 <7 /^L/V^, b° h ^180; (10) 

and for <r(mhos/m)// M H»<<(€/2) X10" 4 : 



K v o-1.7778XlO- 2 C el^-if]& nf 6^90, (11) 
2^1.7778X10- f ayV^/i&, K^90. (12) 
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Figure 2. The parameter b* or b° h . 
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Figure 3. The parameter K v ' or K h . For general A*, 
K(k) = CoK(4:l3). 



Curves of b° and K versus / M hz are shown in 
figures 2 and 3 for various combinations of e and a 
corresponding to poor ground, average ground, good 
ground, and sea water. A value of &=4/3 was 
assumed in obtaining the curves of figure 3 ; however, 
they may be used for other k through the simple 
relationship: K(k) = C Q K(4/3). 

The parameter /3 appearing in (2) is plotted as a 
function of K and b° in figure 4. The limiting value 
of |8 =1.607 for if^O.Ol may be used for almost 
all practical cases of horizontal polarization. On 



the other hand for vertical polarization /3 may 
range from 1.607, corresponding to high frequency 
propagation, down to 0.700, corresponding to low 
frequency propagation over sea water. (See fgs. 
3 and 4.) 

With K, b°, and O evaluated, the groundwave 
attenuation may now be calculated by (1) and 
figures 5 and 7. C X (K, b°), obtained from figure 5, 
has the limiting values <7 X (0, 6°) = 20.03 for K->0 
and Ci(a>, b°) = 20.94 for K-^<» ; for any value of K 
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Figure 4. The parameter /3 . 



it can be seen that C x does not vary more than about 
6dB. 

The functions O(x ), F(x{), and F(x 2 ) are shown in 
figure 7, G(x ) being the uppermost curve and defined 
by: 

G(xo) =5.751X10- 2 a:o-10 log x . (13) 

In general the height function F{x) depends on K 
and 6°, except for large values of x L or x 2 , in which 



case F(x)^G(x). For Xi or x 2 approaching zero, the 
limiting values of F(x) are plotted in figure 6 which 
may thus be used to help interpolate for those values 
of K and 6° not shown in figure 7. 

Since (1) is based upon only the first term of the 
residue series, a fairly good indication of its range of 
validity may be gained from the ratio of the second 
term T 2 to the first term T x of the series. If we 
require the error in A t to be less than some given 
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Figure 6. The height gain function F(x=0) for zero antenna heights. 



value, say 8, we have that: 

20log|l + T 2 /T 1 |<20log[l + |T 2 /T 1 |]<5. 



(14) 



This is the procedure used to obtain the A(x 1(2 ) 
curves shown in the lower right hand corner of 
figure 7. Using the limiting values of p (1.607 and 
0.700), then for A t to be accurate to within 1.5 dB 
(approximately) it may be shown that: 

x -x 1 A(xi)-x 2 A(x 2 )>33d f (00=1.607, K< 0.1) (15a) 

XQ-z 1 A(x 1 )-X2A(x 2 )>n5, (/3o=0.7O0, #>10). (15b) 

Notice that in certain cases (1) provides a good 
approximation to the attenuation even for just 
within line-of-sight paths as long as the appropriate 
condition (15a) or (15b) holds. For values of O 
lying between the two limits, linear interpolation 
between the A(x) curves of figure 7 and the two 
minimum values of (15a) and (15b) will give a fair 
approximation to the range of validity of (1). 



Using linear interpolation: 

X — XiA(x h $o) — X 2 A(x 2) ft))>£mJn, ( 16a ) 

where 



2w n =335-242.6(1.607-A)) 
A(x, A>) = A(z, 1.607) + 1.103(1 
X{A(x, 0.700)-A 



i(1.607 -&) V 
(x, 1.607)}; J 



(16b) 



A(x, 0.700) and A(x, 1.607) are the values read from 
the upper and lower curves respectively in figure 7. 

3. Examples 

In order to use (1) to calculate the field strength E 
in a particular situation* a reference field E must be 
specified. The relationship between A t and E is 
then simply: 



E=E -A t . 



(17) 
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For instance, the reference field (in dB above 1/xV/m) 
assumed for the curves presented in the CCIR 
Atlas [Atlas, 1955] is 

E =20 log [2.22X10 5 /rf (km)] 

= 106.93-20 log rfo(km). (18) 

The following examples are given to illustrate the 
use of (1) in calculating field strength. 

Given a particular propagation path in which 
&=4/3, e = 80, (7=4 mhos/m, /mhz = 30, d =W0 km, 
A,i = 100 m, and A, 2 =1000 m, then from (5), di=41.225 
km and d 2 = 130.365 km; this corresponds to a within 
line-of-sight path. The fields for both vertical and 
horizontal polarization are then calculated as follows: 

(1) Vertical polarization. From (9) the param- 
eters K v and 6° are calculated as K v ^0.28, b° v ^0, 
and from the curves of figures 4 and 5, /5 = 1.38, 
Ci(0.28, 0) = 19.2. From (2) the x's are calculated 
as £o=686, x\ = lll , x 2 = 559, and the criterion given 
by (16) indicates that (1) may be used in the present 
situation: 

686-177(0. 16)-559(0.62) 

= 311>280, from (16a) and (16b). 

Then from (1) and figure 7, 

^=11.0+25.5+0.5-19.2 = 17.8. 

To compare with the CCIR Atlas, (17) and (18) are 
used giving the field strength as 

£=62.8-17.8=:45.0 (dB above 1 juV/m). 

The value as read from the Atlas is £'=44.4. 

(2) Horizontal polarization. From (10) and fig- 
ures 4 and 5, K h ^12Xl0~\ 6°^ 180, &= 1.607, 
<7i(l(T 4 , 180)=20.0; from (2) z =799, ^=206, 
:r 2 = 651, and from (1) and figure 7, 

^=17.0+23.5—4.5-20.0=16.0. 

Equations (17) and (18) give #=62.8 — 16.0=46.8, 
whereas the value as read from the Atlas is £'=46.0. 
Radio propagation studies are often presented in 
terms of the propagation loss L v of a radio system 
[Norton, 1959; Wait, 1959]. The relationship be- 
tween attenuation and propagation loss is 

i,=32.45+20log/ M Hzrfo(km)+^-^, (19) 

where 6 P denotes the path antenna gain which, for 
negligible polarization coupling loss, is just the sum 



of the free space antenna power gains above an 
isotropic antenna. 

The system loss L s of a radio communication 
system is defined as the decibel ratio of the power 
input to the terminals of the transmitting antenna 
P t and the resultant power available at the terminals 
of the receiving antenna P a : 



L s =P t -P a . 



(20) 



Propagation loss and the practicably measurable 
system loss are then related by 



J^s — P t Pj. — lj v -\-L t -\-Li r , 



(21) 



where L t (transmitter) and L r (receiver) are the 
decibel ratios of the actual antenna input resistances 
r to their free space radiation resistances r f : 



L t>r =10 log (r/r f ). 



(22) 



Curves of L t and L r for electric dipoles and small 
loop antennas without ground screens over finitely 
conducting ground are available in the literature 
[Vogler and Noble, 1963]. The effects of ground 
screens are considered in papers by Wait and Surtees 
[1954] and by Wait [1956]. 
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